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ORIGINAL ARTICLE 

Network inefficiencies in autism spectrum disorder at 
24 months 

JD Lewis\ AC Evans\ JR Pruett^ K Botteron^ L Zwaigenbaum^ A Estes^ G Gerig^ L Collins\ P Kostopoulos\ R McKinstry^ S Dager^ 
S Paterson^ RT Schultz^ M Styner^'^ H Hazlett^ and J Piven^ for the IBIS network 

Autism spectrum disorder (ASD) is a developmental disorder defined by behavioral symptoms that emerge during the first years of 
life. Associated with these symptoms are differences in the structure of a wide array of brain regions, and in the connectivity 
between these regions. However, the use of cohorts with large age variability and participants past the generally recognized age of 
onset of the defining behaviors means that many of the reported abnormalities may be a result of cascade effects of 
developmentally earlier deviations. This study assessed differences in connectivity in ASD at the age at which the defining 
behaviors first become clear. There were 1 1 3 24-month-old participants at high risk for ASD, 31 of whom were classified as ASD, and 
23 typically developing 24-month-old participants at low risk for ASD. Utilizing diffusion data to obtain measures of the length and 
strength of connections between anatomical regions, we performed an analysis of network efficiency. Our results showed 
significantly decreased local and global efficiency over temporal, parietal and occipital lobes in high-risk infants classified as ASD, 
relative to both low- and high-risk infants not classified as ASD. The frontal lobes showed only a reduction in global efficiency in 
Broca's area. In addition, these same regions showed an inverse relation between efficiency and symptom severity across the high- 
risk infants. The results suggest delay or deficits in infants with ASD in the optimization of both local and global aspects of network 
structure in regions involved in processing auditory and visual stimuli, language and nonlinguistic social stimuli. 
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INTRODUCTION 

Autism spectrum disorder (ASD) is a developmental disorder 
defined by impairments in social communication and social 
interaction and a restricted repertoire of activities and interests.^ 
Recent research has indicated that abnormalities in connectivity 
may underlie many of the behavioral symptoms. This research, 
however, is equivocal with respect to both the nature of 
connectivity abnormalities and the regions in which they occur, 
and currently provides little insight into their early developmental 
origins. Numerous studies have reported long-range under- 
connectivity in ASD,^"^ some have reported short-range over- 
connectivity accompanying long-range under-connectivity,^"^^ 
and some have reported a less clearly divided mixture of over- 
and under-connectivity.^^"^^ Such connectivity abnormalities have 
been variously reported for virtually all pairs of brain 
regions.^-^'^^-^° 

This disparate set of findings reflects a variety of important 
methodological differences between studies, as well as small 
sample sizes and substantial etiologic and clinical heterogeneity. 
Prominent among the sources of variability is the wide variety of 
ages of the participants in these various studies. ASD is a 
developmental disorder involving complex networks comprising 
multiple spatially distributed brain regions which are developing 
in the context of a clearly altered growth trajectory. Children with 
ASD show abnormally rapid brain growth during the first years of 



life,^^"^^ and abnormalities in the development of connectivity.^^ 
These brain-based abnormalities are likely to contribute to a 
developmental cascade of behavioral and neural abnormalities.^^ 
Thus, differences in findings of connectivity abnormalities in ASD 
may be owing to differences in the ages of the participants in 
various studies. The multitude of abnormalities in connectivity 
reported in adults with ASD, or even children with ASD, may tell us 
little about where in the brain these abnormalities first emerge. 
Identifying where connectivity abnormalities originate and how 
they develop has the potential to inform as to their cause, and to 
suggest more etiologically relevant paths to successful interven- 
tion. However, if we are to avoid confusing downstream 
developmental effects with source effects, this identification must 
be on the basis of data from individuals with ASD as early in 
development as possible, and with as little variation in age as 
possible. 

Behavioral abnormalities may appear during the first year of life 
or early in the second year, but typically the characteristic 
symptoms of ASD necessary for classification only become stable 
by the end of the second year.^^"^^ If connectivity abnormalities 
underlie these symptoms of ASD, then differences in connectivity 
should also be present by the end of the second year of life, but 
downstream developmental effects should be minimal. Thus, to 
localize the early neural abnormalities in ASD, we assess 
connectivity at 24 months of age. 
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We utilize a whole-brain network analysis to assess white matter 
connectivity abnormalities in three groups of 24-month-olds: 
those with symptoms of ASD who also have an older sibling with 
ASD; those having an older sibling with ASD, but who are not 
themselves developing symptoms of ASD; and typically develop- 
ing 24-month-olds with an older sibling and no first-degree 
relative with ASD or intellectual disability. Children with a sibling 
who has been diagnosed with ASD, but who do not receive a 
diagnosis themselves, often exhibit milder ASD symptoms.^^'^^ 
Thus this group might be expected to show a milder form of any 
connectivity abnormalities found in infants with significant 
symptoms of ASD. 

A network analysis assesses connectivity on a whole-brain basis, 
rather than in terms of individual connections, thus providing 
measures which are relevant to overall brain organization.^^ Brain 
areas are treated as nodes in a graph; the connections between 
the brain areas constitute the edges of the graph. Properties of the 
overall network, or of subnetworks, associated with each 
participant can then be measured, and group differences 
assessed. This analytic approach has proved effective in char- 
acterizing differences in brain networks in a variety of clinical 
populations, including ASD.^'^^ 

Network analysis methods have evolved over the past 15 years, 
from straightforward applications of graph theory, which assess 
only network topology,^^ to more sophisticated approaches which 
take into account the spatial aspects of connectivity to assess the 
efficiency of information transfer within the network.^^"^^ Such 
approaches utilize measures of the length and strength of 
connections between all pairs of anatomical regions to estimate 
how efficiently information can be transferred between regions. 
This provides a way to assess the impact of abnormalities in spatial 
organization on communication, both across the brain and within 
local subnetworks. We use probabilistic tractography to estimate 
the strength of connectivity between all pairs of regions, and the 
length of the connections between regions. We then assess the 
group differences in local and global efficiency for each brain 
region, as well as the relation of efficiency to symptom severity, to 
determine where communication efficiency might be impaired in 
the brains of 24-month-olds with symptoms of ASD. 

MATERIALS AND METHODS 

Participants 

The participants for this study were drawn from an ongoing multisite study 
of brain and behavioral development in infants at high risk for autism due 
to having an older sibling with ASD, known as the IBIS (Infant Brain 
Imaging Study), a study funded by the National Institute of Health via an 
Autism Center of Excellence grant. IBIS has enrolled and collected 
neuroimaging and behavioral data from these high-risk infants at 6, 12 
and 24 months of age, with confirmation upon enrollment of the diagnosis 



of their older sibling by medical record review and the Autism Diagnostic 
Interview. Data are also collected for infants deemed at relatively low risk 
for autism by virtue of having no first-degree relative with ASD or 
intellectual disability (as assessed by the Social Communication Ques- 
tionnaire and the Family Interview for Genetic Studies), and an older 
sibling. A parent or legal guardian gave informed consent. Exclusionary 
criteria include the following conditions: (1) diagnosis or physical signs of 
genetic conditions/syndromes (for example, fragile X syndrome), (2) 
significant medical or neurological conditions affecting growth, develop- 
ment or cognition (for example, CNS infection, seizure disorder, congenital 
heart disease), (3) sensory impairments such as vision or hearing loss, (4) 
low birth weight ( < 2000) or prematurity ( < 36 weeks gestation), (5) 
evidence of significant perinatal brain injury or exposure to in-utero 
neurotoxins (for example, alcohol, selected prescription medications), (6) 
non-English speaking families, (7) contraindication for MRI (magnetic 
resonance imaging; for example, metal implants), (8) adopted children and 
(9) family history of first-degree relative with intellectual disability (not 
associated with ASD), psychosis, schizophrenia or bipolar disorder. The 
data were collected at four sites (University of North Carolina at Chapel Hill, 
University of Washington, Children's Hospital of Philadelphia, Washington 
University in St Louis); quality control and pre-processing were done at the 
University of Utah, the University of North Carolina at Chapel Hill and the 
Montreal Neurological Institute (McGill University), the Data Coordinating 
Center. The study was approved by the institutional review boards at 
all sites. 

The subset of data used in the current study includes all 24-month data 
that passed quality control as of June 201 2. These data divide according to 
whether the infants are high-risk or low-risk, and according to whether 
they showed ASD symptoms, as indexed by the Autism Diagnostic 
Observation Schedule (ADOS),'^° above the clinical cutoff for ASD. 
Participants were administered the ADOS at 24 months of age by trained 
researchers. Descriptives for the participants are shown in Table 1, with 
ADOS scores transformed to calibrated severity scores, as per Gotham 
et al.^^ Note that Table 1 includes high-risk infants above the clinical cutoff 
for ASD (HR^°^), high-risk infants below the clinical cutoff for ASD (HR^^^) 
and low-risk infants below the clinical cutoff for ASD (LR"^^*^). There were 
three low-risk infants who received ADOS scores above the cutoff for ASD; 
their data were excluded. All other infants (a total of 136) were included in 
the analysis. 

Imaging and image processing 

MRI scans were done at each site on Siemens 3T TIM Trio scanners 
(Siemens Medical Systems, ErIangen, Germany) with 12-channel head coils. 
Intra- and inter-site reliability was verified.^^ Data were collected while 
infants were naturally sleeping. Three types of images were acquired: (1) 
Trweighted images using a 3D MPRAGE sequence 
(resolution = 1.0 X 1.0 X 1.0 mm; TE = 3.16ms; TR = 2400 ms; matrix = 224 x 
256); (2) r2-weighted images using a 3D FSE sequence 
(resolution = 1.0x1.0x1.0 mm; TE = 499 ms; TR = 3200 ms; matrix = 256 x 
256); and (3) diffusion-weighted images using a 2D echo planar sequence 
(resolution = 2.0x2.0x2.0; TE= 102ms; TR= 12 800ms; directions = 25; 
b-values: unique values evenly distributed between 0 and lOOOsmm"^). 
Data were quality controlled, starting at the MRI console; where motion 
artifacts were observed, the scan sequence was repeated to possibly 
acquire artifact free data. 

The Tr and r2-weighted images were subjected to manual quality 
control, and corrected for geometric distortion.^^ The diffusion-weighted 
images were cleaned of motion and other artifacts using DTIPrep,^^ which 
corrects artifacts where possible, and excludes directions from the data 
when correction is not possible. Only data sets with acceptable 7"r, T2- and 
diffusion-weighted images were included in the analysis. 

The distortion-corrected 7"r and r2-volumes were then processed with 
CIVET, a fully automated structural image analysis pipeline developed at 
the Montreal Neurological Institute. CIVET corrects intensity nonuniformi- 
ties using N3;^^ aligns the input volumes to the Talairach-like ICBM-152-nl 
template;"^^ classifies the image into white matter, gray matter, cere- 
brospinal fluid and background;"^^"^^ extracts the white matter and pial 
surfaces;^^ and warps these to a common surface template.^° The CIVET 
results were used to construct the seed, stop and target masks for use with 
FSL's probtrackx.^^ Seed masks specify the voxels from which tracts 
propagate; seed masks were the entire white matter. Stop masks specify 
where tract propagation is halted; stop masks were the voxels on the 
boundary of the white matter. Target masks determine the interpretation 
of tracts; target masks were the cortical labels of the automatic anatomical 



Table 1. Sample characteristics 


Age (months) 


















N 


Mean s.d. 


N 


Mean 


s.d. 


N 


Mean s.d. 


Male 
Female 


24 
7 


24.85 0.64 
25.05 0.82 


51 
31 


24.76 
24.72 


0.74 
1.11 


12 
11 


24.62 0.70 
24.45 0.54 


ADOS severity 
Male 
Female 




6.13 1.60 

5.14 1.86 




1.36 
1.29 


0.70 
0.64 




1.33 0.65 
1.22 0.67 


Abbreviations: HR^°^ high-risk infants with ADOS scores above the clinical 
cutoff for ASD; HR"^^^, high-risk infants with ADOS scores below the clinical 
cutoff for ASD; LR"^^^, low-risk infants with ADOS scores below the clinical 
cutoff for ASD. 
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labelling (AAL) atlas,^^ a parcellation based on anatomical landmarks, 
predominately sulci, defined on the common surface template.^^ 

The diffusion-weighted images were unwarped via nonlinear registra- 
tion to the distortion-corrected r2-volume, and then preprocessed for 
probabilistic tractography with FSL's bedpostx.^^ The default parameters 
for bedpostx were used: a maximum of two fibers per voxel, a 
multiplicative factor of 1 on the prior for the second fiber and a burn-in 
of 1000. The unwarped diffusion volumes were also affine registered to the 
Trvolumes using FSL's flirt^"^ to provide the mapping from the seed, stop 
and target masks to the diffusion space. Probabilistic tractography, utilizing 
FSL's probtrackx was then seeded twice from 10 000 random locations 
within each voxel of the seed masks: with distance-bias correction^^'^^ to 
generate the number of tracts connecting voxels in the target mask, and 
without distance-bias correction to generate the physical lengths of those 
tracts. These results were then compiled for each AAL region, generating 
undirected matrices of the total number of connections between each pair 
of AAL regions, and the mean physical length of those connections. The 
total number of connections between each pair of AAL regions divided by 
the average surface area of the two AAL regions is referred to as 
connection 'strength'. 

Analysis 

Network analyses investigate the organization of sets of connections, 
rather than considering the connections independently. The metric utilized 



here is efficiency, as defined by Latora and Marchiori." 



Network 



efficiency is the capacity to exchange information across a network. Latora 
and Marchiori^^ defined the efficiency £,y in the communication between 
nodes / and j to be inversely proportional to the shortest path length d,y 
between nodes / and j. They take the shortest path length c/,y to be the 
smallest sum of the physical distances throughout all of the possible paths 
from / to j in the graph. The efficiency of a network, G, is then 

F(G) = ^'^j^^^'j = 'I J_ 
^ ^ N{N-^) N{N -^)^i^j^(^clij 

where N is the number of nodes in the network graph G; £,y is the efficiency 
of the connection between nodes / and j; and d,y is the length of the 
shortest path, in terms of physical distances, between nodes / and j. This 
measure is normalized by E(Gideal)> the fully connected network. This 
formulation is importantly different from most graph theoretic measures, 
which assess only the topological properties of networks;^"^'^^'^^ the 
measures of efficiency take into account the physical distances involved in 
information transfer, and so presumably relate more closely to the 
neurobiological substrates. 

Latora and Marchiori^^ apply this formulation to both the entire network, 
which they refer to as global efficiency, and to the subnetworks of the 
immediate neighbors of each node; they define local efficiency as the 
mean of E(G,), for all nodes /, where G, is the subgraph of all the neighbors 
of node /. These definitions give a single measure of local efficiency and of 
global efficiency for the entire network. But, the definitions can be 
reformulated to provide measures of efficiency for each node, or for 
collections of nodes. Achard and Bullmore^^ define nodal efficiency, which 
we will refer to as nodal global efficiency, as the inverse of the harmonic 
mean of the minimum path length between a node, /, and all other nodes 
in the network. Utilizing the physical distances, as per Latora and 
Marchiori,^^ the nodal global efficiency of node / is thus 



fnodalglobal(<J, ■ 



(A/ - 1) {N -^)^j^GJ^jcli 



—y 



1 



where N is the number of nodes in the network graph G; £,y is the efficiency 
of the connection between nodes / and j; and d,y is the length of the 
shortest path, in terms of tractography-based measures of physical 
distances, between nodes / and j. The definition of local efficiency can 
likewise be parsed to provide a measure of nodal local efficiency; recall 
that the local efficiency of a network is the mean of E(G,), for all nodes /, 
where G, is the subgraph of all the neighbors of node /. Thus the nodal 
local efficiency of node / is simply 



£'nodallocal(G, /) -- 



1 



A/g, (A/g,. - 1 ) A/g,. (A/g,. - 1 )^j^k e G;clj 



1 



fy/c 



where A/g, is the number of nodes in the subgraph G, consisting of all of the 
neighbors of/; Sjk is the efficiency of the connection between nodes j and k; 
and djk is the length of the shortest path, in terms of tractography-based 
measures of physical distances, between nodes j and k. 



These definitions treat connections in a binarized fashion. But the 
strengths of the connections in the brain are critical to an accurate 
assessment of its efficiency.^^ We thus utilize a version of these measures 
that does not discard information about connection strength by thresh- 
olding, but rather incorporates connection strengths. On the basis of 
Rubinov and Sporns,^^ we define the weighted distance between nodes / 
and j as 

where c/,y is the length of the shortest path, in terms of tractography-based 
measures of physical distances, between nodes / and y; and Wjj is the 
connection strength between nodes / and j. The length of the shortest path 
in terms of weighted distances is thus related to the travel time along the 
path, decreasing with the strengths of the connections involved, and 
increasing with the tractography-based measures of the physical distances 
involved. Our weighted formulations of nodal global efficiency and nodal 
local efficiency, also based on Rubinov and Sporns,^^ are 



r-weighted 
nodalglobal 



1 



(N-1)^ 

where N is the number of nodes in the network graph G; andc/jj^ is the 
shortest path, in terms of weighted distance, between nodes / and j; and 



r-weighted 



(G,0 = 



1 

A/g,(A/g,-1) 



1/3 



where A/g, is the number of nodes in the subgraph G, consisting of all of the 
neighbors of /; c/])J is the shortest path, in terms of weighted distance, 
between nodes j and k; and w,y and Wj^ are the connection strengths 
between nodes / and j, and / and k, respectively. As per Latora and 
Marchiori,^^ these measures are normalized by considering the fully 
connected network. 

Group differences in nodal local efficiency and nodal global efficiency 
were assessed via statistical linear models, controlling for age, sex and site, 
in models which included all three groups. To determine the extent to 
which any group differences in network structure were driven by 
individual behavioral differences, we also assessed the relation between 
the measures of network efficiency and the ADOS-based calibrated 
severity scores^^ within the high-risk infants, controlling for age, sex and 
site. To gain further insight into the proper interpretation of any group 
differences in nodal local efficiency and nodal global efficiency, a post hoc 
analysis was performed to assess group differences in three additional 
properties of the networks associated with each node: the strengths of the 
connections comprising the paths contributing to the efficiency measures; 
the tract-based lengths of those paths; and the number of edges in those 
paths. For each of these measures, the contribution associated with the 
path between nodes j and k to the nodal local measure at node / is 
determined by the connection strength between nodes / and j, and / and k, 
exactly as in the case of nodal local efficiency. The nodal global measures 
at node / are the mean of the measures associated with the paths 
contributing to nodal global efficiency at node /. 



RESULTS 

Relative to LR"^^^ infants, the HR''^^ infants showed reductions in 
nodal local efficiency bilaterally in the temporal and occipital lobes 
in inferior and medial regions, and predominately in the right 
hemisphere for lateral regions, extending to the supramarginal 
gyrus (Figure 1, top left). Relative to LR"^^^ infants, the HR^°^ 
infants showed reductions in nodal global efficiency in bilateral 
temporal and occipital regions, extending in the right hemisphere 
to the angular and supramarginal gyri; the left hemisphere 
reductions include the pars triangularis and orbital gyrus. 
(Figure 1, bottom left). There were no significant increases in 
either nodal local or nodal global efficiency in the HR'^°^ infants 
relative to the LR"^^^ infants. The group differences in nodal local 
efficiency for HR'^^^ compared with HR"^^^ infants are similar to 
those in the comparison of the V\\f^^ and LR"^^^ infants, but left 
lateralized and less extensive. Relative to the HR"^^*^ infants, the 
HR'^°^ infants showed reductions in nodal local efficiency over the 
left hemisphere temporal and occipital lobes, and the precuneus; 
the right hemisphere differences are mostly nonsignificant 
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Figure 1. The connparison of nodal local and nodal global efficiency for LR"^^*^ versus HR^°^ infants, and for HR'^^*^ versus HR'^°^ infants, in a 
nnodel with all three groups, controlling for age, sex and site. The top half of the figure shows the f- and p-statistic for the connparison of nodal 
local efficiency in each region of the AAL atlas; the bottonn half of the figure shows the t- and p-statistic for the connparison of nodal global 
efficiency. The left colunnn shows the connparisons for LR"^^*^ versus HR^°^ infants; the right colunnn for HR'^^*^ versus HR infants. A positive f- 
statistic represents a relative reduction in efficiency in HR^°^ infants. No regions show significantly increased nodal local or global efficiency in 
HR^°^ infants. Significant reductions in nodal local efficiency in HR^°^ connpared with LR"^^*^ infants are seen bilaterally in the tennporal and 
occipital lobes in inferior and nnedial regions, and predonninately in the right hennisphere for lateral regions, extending to the suprannarginal 
gyrus. Significant reductions in nodal local efficiency in HR^°^ connpared with HR'^^ infants are seen in the left tennporal and occipital lobes, 
extending into the precuneus, and in several posterior regions in the right hennisphere. Significant reductions in nodal global efficiency in 
HR^°^ infants connpared with LR"^^*^ infants are seen bilaterally in the tennporal and occipital lobes, extending in the right hennisphere to the 
angular and suprannarginal gyri. The left pars triangularis and nnedial orbital gyrus also show reduced nodal global efficiency. Significant 
reductions in nodal global efficiency in HR infants connpared with HR'^^*^ infants are seen in bilateral tennporal lobes, extending in the left 
hennisphere to the occipital cortex, and in the right hennisphere to the angular and suprannarginal gyri. The left hennisphere reductions include 
the nnedial tennporal and occipital regions, and the pars triangularis. The right hennisphere reductions include the precuneus. Note that in all 
cases the reductions are predonninately posterior. HR'^^*^, high-risk infants with ADOS scores below the clinical cutoff for ASD; HR^°^ high-risk 
infants with ADOS scores above the clinical cutoff for ASD; LR"^^*^, low-risk infants with ADOS scores below the clinical cutoff for ASD. 
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(Figure 1, top right). The group differences in nodal global 
efficiency for HR'^^^ compared with HR"^^^ infants are also similar 
to those in the comparison of the HR^°^ and LR"^^^ infants, but 
with greater reductions in the left hemisphere. Significant 
reductions in HR'^^^ infants compared with HR"^^*^ infants are 
seen in bilateral temporal lobes, extending in the left hemisphere 
to the occipital cortex, and in the right hemisphere to the angular 
and supramarginal gyri. The left hemisphere reductions include 
the pars triangularis. The right hemisphere reductions include the 
precuneus (Figure 1, bottom right). There were no significant 
increases in either nodal local or global efficiency in the HR'^^^ 
infants relative to the HR^^^ infants. The HR^^^ and LR^^^ infants 
showed no group differences; but note that the results are 
spatially graded: HR'^^^ infants showed more extensive reductions 
relative to LR"^^^ infants than to the HR"^^^ infants. Also note that 
in all the cases the reductions are predominately posterior. 

The relationship between the measures of network efficiency 
and the ADOS-based calibrated severity scores within the high-risk 
infants are similar to the group differences in efficiency for the two 
subgroups of the high-risk infants, that is, the HR compared 
with HR"^^^ infants. There is a significant inverse relation of severity 
to nodal local efficiency in the left temporal and occipital lobes 
(Figure 2, top). This inverse relation is significant in the lingual 
gyrus, and throughout the temporal lobe, with the exception of 
the middle temporal gyrus. There is a significant inverse relation of 
severity to nodal global efficiency in the left temporal lobe and in 
Broca's area (Figure 2, bottom). This inverse relation is significant 
throughout the temporal lobe, with the exception of the inferior 
temporal gyrus. Note that where there is a significant relation 
between efficiency and severity, in all cases this is an inverse 
relation, thus greater symptom severity corresponds to lesser 
efficiency. Also note that, as with the group differences, this 
inverse relation is only present in posterior regions, and in Broca's 
area for nodal global efficiency. 

The post hoc analyses of the group differences in the strengths 
of the connections comprising the paths contributing to the 
efficiency measures, the tract-based lengths of those paths, and 
the number of edges in those paths indicates that these results 
stem from both weaker connections and abnormalities in network 
structure. The reductions in nodal local and nodal global efficiency 
in HR'^°^ infants are partially paralleled by reductions in 
connection strength in the paths involved. Relative to LR"^^^ 
infants, the HR^^^ infants showed reductions in nodal local 
connection strength bilaterally in the temporal, parietal and 
occipital lobes, and in the left medial orbital gryus (Figure 3, top 
left). Relative to LR"^^^ infants, the HR^°^ infants showed 
reductions in nodal global connection strength in the right 
inferior temporal and occipital lobes (Figure 3, bottom left). There 



Figure 2. The relation of ADOS-based calibrated severity to network 
efficiency in the high-risk infants, controlling for age, sex and site. 
The top half of the figure shows the f- and p-statistic for the relation 
of severity and nodal local efficiency in each region of the AAL atlas; 
the bottom half of the figure shows the t- and p-statistic for the 
relation of severity and nodal global efficiency. A negative f-statistic 
represents an inverse relation between severity and efficiency, that 
is, greater symptom severity corresponds to lesser efficiency. There 
is a significant inverse relation of severity to nodal local efficiency in 
the left lingual gyrus, and throughout the left temporal lobe, with 
the exception of the middle temporal gyrus. There is a significant 
inverse relation of severity to nodal global efficiency in Broca's area, 
and throughout the left temporal lobe, with the exception of the 
inferior temporal gyrus. Note that all significant relations are inverse 
relations, and that these are only present in the left temporal and 
occipital lobes, and in Broca's area for nodal global efficiency. AAL, 
automatic anatomical labelling; ADOS, Autism Diagnostic Observa- 
tion Schedule. 



were no significant increases in either nodal local or nodal global 
connection strength in the HR^°^ infants relative to the LR"^^^ 
infants. Relative to HR"^^^ infants, the HR^°^ infants showed 
reductions in nodal local connection strength in bilateral temporal 
and parietal lobes, the right precuneus and gyrus rectus, and the 
left pars triangularis and medial orbital gyrus (Figure 3, top right). 
Relative to HR^^^ infants, the HR''^^ infants showed reductions in 



Efficiency and Severity 





© 2014 Macmillan Publishers Limited 



Translational Psychiatry (2014), 1-11 




0.0 0.05 0.0 0.05 



Figure 3. The comparison of the mean connection strength in the paths comprising the measures of nodal local and nodal global efficiency 
for LR"^^*^ versus HR^°^ infants, and for HR'^^*^ versus HR infants, in a model with all three groups, controlling for age, sex and site. The top 
half of the figure shows the t- and p-statistic for the comparison of nodal local connection strength in each region of the AAL atlas; the bottom 
half of the figure shows the f- and p-statistic for the comparison of nodal global connection strength. The left column shows the comparisons 
for LR"^^*^ versus HR^°^ infants; the right column for HR"^^*^ versus HR infants. A positive f-statistic represents a relative reduction in 
connection strength in HR^°^ infants. No regions show significantly increased connection strength in HR^°^ infants. Significant reductions in 
nodal local connection strength in HR'^°^ compared with LR"^^*^ infants are seen bilaterally in the temporal, parietal and occipital lobes, and in 
the left medial orbital gryus. Significant reductions in nodal local connection strength in HR^°^ compared with HR"^^*^ infants are seen in 
bilateral temporal and parietal lobes, the right precuneus and gyrus rectus, and the left pars triangularis and medial orbital gyrus. Significant 
reductions in nodal global connection strength in HR^°^ infants compared with LR"^^*^ infants are seen in the right fusiform gyrus and inferior 
occipital lobe. Significant reductions in nodal global connection strength in HR^°^ infants compared with HR"^ infants are seen in the right 
temporal and occipital lobes. Note that, as with efficiency, in all cases the reductions are predominately posterior. HR'^^*^, high-risk infants with 
ADOS scores below the clinical cutoff for ASD; HR^°^ high-risk infants with ADOS scores above the clinical cutoff for ASD; LR"^^^, low-risk 
infants with ADOS scores below the clinical cutoff for ASD. 
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Figure 4. The connparison of the nnean length of the paths connprising the nneasures of nodal local and nodal global efficiency for LR"^^*^ versus 
HR^°^ infants, and for HR'^^*^ versus HR^°^ infants, in a nnodel with all three groups, controlling for age, sex and site. The top half of the figure 
shows the t- and p-statistic for the connparison of nodal local nnean path length in each region of the AAL atlas; the bottonn half of the figure 
shows the f- and p-statistic for the connparison of nodal global nnean path length. The left colunnn shows the connparisons for LR"^^*^ versus 
HR^°^ infants; the right colunnn for HR versus HR'^°^ infants. A positive f-statistic represents a relative reduction in nnean path length in 
HR^°^ infants. Significant reductions in nodal local nnean path length in HR^°^ connpared with LR"^^*^ infants are seen bilaterally in the 
tennporal, parietal and occipital lobes, and in the left nnedial orbital gryus. Significant reductions in nodal local nnean path length in HR^°^ 
connpared with HR'^^*^ infants are seen in bilateral tennporal, parietal and occipital lobes, the right precuneus, paracentral lobule and gyrus 
rectus, and the left pars triangularis and nnedial orbital gyrus. No regions show significantly increased nodal local nnean path length in HR'^°^ 
infants relative to either LR"^ or HR"^^*^ infants. Significant increases in nodal global nnean path length in HR^°^ infants connpared with LR"^^*^ 
infants are seen in the right parietal lobe, the right pars triangularis and the left medial orbital gyrus. Significant increases in nodal global 
nnean path length in HR infants connpared with HR *^ infants are seen bilaterally in the parietal lobe, and in the right tennporal lobe and left 
nnedial orbital gyrus. No regions show significantly decreased nodal global nnean path length in HR^°^ infants relative to either LR"^^*^ or HR"^^*^ 
infants. Note that, as with efficiency, in all cases the group differences are predonninately posterior. AAL, autonnatic anatonnical labelling; 
HR"^^*^, high-risk infants with ADOS scores below the clinical cutoff for ASD; HR^°^ high-risk infants with ADOS scores above the clinical cutoff 
for ASD; LR"^^*^, low-risk infants with ADOS scores below the clinical cutoff for ASD. 
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Figure 5. The comparison of the mean number of edges in the paths comprising the measures of nodal local and nodal global efficiency for 
LR"^^*^ versus HR infants, and for HR'^^*^ versus HR infants, in a model with all three groups, controlling for age, sex and site. The top half 
of the figure shows the t- and p-statistic for the comparison of nodal local mean number of edges in each region of the AAL atlas; the bottom 
half of the figure shows the t- and p-statistic for the comparison of nodal global mean number of edges. The left column shows the 
comparisons for LR"^^*^ versus HR'^°^ infants; the right column for HR'^^*^ versus HR^°^ infants. A positive f-statistic represents a relative 
reduction in mean number of edges in HR^°^ infants. Significant reductions in nodal local mean number of edges in HR compared with 
LR"^^*^ infants are seen bilaterally in the temporal, parietal and occipital lobes, and in the left medial orbital gryus. Significant reductions in 
nodal local mean number of edges in HR compared with HR infants are seen in bilateral temporal and parietal lobes, the right 
precuneus, paracentral lobule, and gyrus rectus, and the left cuneus, pars triangularis and medial orbital gyrus. No regions show significantly 
increased nodal local mean number of edges in HR^°^ infants relative to either LR"^^*^ or HR"^^*^ infants. Significant increases in nodal global 
mean number of edges in HR'^°^ infants compared with LR"^^*^ infants are seen in bilateral temporal, parietal and occipital lobes, the right 
superior frontal gyrus, and the left medial orbital gyrus. Significant increases in nodal global mean number of edges in HR^°^ infants 
compared with HR"^^*^ infants are seen in bilateral temporal, parietal and occipital lobes, the left superior frontal gyrus, medial orbital gyrus, 
gyrus rectus and pars opercularis, and the right pars triangularis, paracentral lobule and cinqulate cortex. No regions show significantly 
decreased nodal global mean number of edges in HR^°^ infants relative to either LR"^^*^ or HR infants. Note that, as with efficiency, in all 
cases the group differences are predominately posterior. HR"^^*^, high-risk infants with ADOS scores below the clinical cutoff for ASD; HR^°^ 
high-risk infants with ADOS scores above the clinical cutoff for ASD; LR"^^*^, low-risk infants with ADOS scores below the clinical cutoff for ASD. 
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nodal global connection strength in the right temporal and 
occipital lobes (Figure 3, bottom right). There were no significant 
increases in either nodal local or nodal global connection strength 
in the HR''°^ infants relative to the HR"^^^ infants. Note that, as 
with efficiency, in all cases the reductions are predominately 
posterior. 

The reductions in nodal local and nodal global efficiency in 
HR'^°^ infants appear to be not only the result of weaker 
connections, but also abnormalities in network structure. Relative 
to LR"^^^ infants, the HR^°^ infants showed significantly reduced 
nodal local mean path length in bilateral temporal, parietal and 
occipital lobes, and the left medial orbital gyrus (Figure 4, top left). 
The fact that this group difference is almost identical to the group 
difference in nodal local mean number of edges (Figure 5, top left) 
indicates that these networks have reduced spatial extent. Relative 
to LR"^^^ infants, the HR'^°^ infants showed significantly increased 
nodal global mean path length in the right parietal lobe, the right 
pars triangularis, and the left medial orbital gyrus (Figure 4, 
bottom left). The comparison of nodal global mean number of 
edges also shows increases in HR'^°^ infants in posterior regions, 
more bilaterally in this case, as well as in the right lateral orbital 
and superior frontal gyri, and the left medial orbital gyrus 
(Figure 5, bottom left). Thus, relative to LR"^^^ infants, HR^°^ 
infants have on average longer, less direct paths connecting a 
node with other brain regions. Relative to HR"^^^ infants, the HR^°^ 
infants showed significantly reduced nodal local mean path length 
in bilateral temporal and parietal lobes, the left occipital lobe, pars 
triangularis and medial orbital gyrus, and the right precuneus, 
paracentral lobule, and gyrus rectus (Figure 4, top right). Again, 
this group difference is exactly mirrored by the group difference in 
nodal local mean number of edges (Figure 5, top right) indicating 
that nodal local networks have reduced spatial extent in HR^^^ 
infants. Relative to HR"^^^ infants, the HR''°^ infants showed 
significantly increased nodal global mean path length in bilateral 
parietal lobes, the right temporal lobe and pars triangularis, and 
the left medial orbital gyrus (Figure 4, bottom right). The 
comparison of nodal global mean number of edges shows 
increases in HR'^^^ infants in bilateral temporal, parietal and 
occipital lobes, as well as in the left medial orbital and superior 
frontal gyri, the left medial orbital gyrus and pars orbitalis, and the 
right pars triangularis, paracentral lobule and cingulate gyrus 
(Figure 5, bottom right). Thus, relative to HR"^^^ infants, HR^°^ 
infants have on average longer, less direct paths connecting a 
node with other brain regions. 

DISCUSSION 

Network efficiency measures from diffusion-based probabilistic 
tractography were significantly reduced, predominately in poster- 
ior brain regions, in HR^°^ infants relative to both LR^^^and HR^^^ 
infants at 24 months of age, and were inversely related to ASD 
symptom severity. HR'^°^ infants showed reduced nodal global 
efficiency, relative to both LR"^^^ and HR"^^^ infants, over the right 
superior temporal, angular and supramarginal gyri, and the left 
pars triangularis, cuneus and lingual gyrus. HR'^^^ infants showed 
reduced nodal local efficiency, relative to both LR^^^ and HR^^^ 
infants, over bilateral inferior temporal gyri, the left fusiform, 
parahippocampal and lingual gyri, temporal pole, and medial 
occipital cortex. Relative to the HR"^^^ infants, the HR''^^ infants 
showed more bilateral reductions in nodal global efficiency and 
more left lateralized reductions in nodal local efficiency. In either 
case, the abnormalities were concentrated in posterior regions; 
the frontal lobes showed no differences, with the exception of the 
anterior portion of Broca's area. The significant differences in the 
HR^^^ infants were all reductions in efficiency. Further, both nodal 
local and nodal global efficiency showed an inverse relation to 
ASD symptoms, as indexed by ADOS-based calibrated severity 
scores within the high-risk infants, over left temporal and occipital 
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regions, and in Broca's area for nodal global efficiency. The regions 
implicated are involved in processing auditory and visual stimuli, 
language and nonlinguistic social stimuli, for example, faces, and 
in memory, all of which show abnormalities in ASD.^^"^^ 

These results align with a number of findings in infant siblings 
of children with autism, and older children with ASD. Infant 
siblings of children with autism have been shown to have 
atypicalities in both visual and auditory processing,^^"^^ as well as 
an inability to integrate audio-visual information/^ in addition, 
connections between the bilateral occipital lobes have been 
implicated in relation to visual orienting in infants who go on to 
develop symptoms of ASD.^^ Studies in toddlers with symptoms of 
ASD have shown reduced correlation in functional data between 
bilateral superior temporal gyri,^^ and reduced activation in 
response to speech stimuli in the left superior temporal gyrus.^^ 
Slightly older children with ASD show reduced correlation in 
functional data between bilateral inferior frontal gyri.^^ Hypoper- 
fusion in bilateral temporal and temporal-parietal regions has also 
been reported in older children with ASD,^^"^^ and a negative 
relation between cerebral blood flow and autism severity in the 
left superior temporal gyrus.^^ Structural abnormalities in older 
children with autism include localized gray matter reductions 
within ventral and superior temporal cortical areas, as well as 
parietal areas.^^ 

Frontal cortico-cortical connectivity abnormalities, which are 
well documented in adolescents and adults with ASD,'''''^'^° 
presumably emerge over development together with the 
behavioral abnormalities with which they are associated. In 
children, the frontal lobes show few network abnormalities, 
whereas posterior regions show widespread abnormalities both 
locally and globally.^'^^ 

The reductions in nodal local and nodal global efficiency shown 
here for the HR'^°^ infants correspond to a number of organiza- 
tional differences. The paths underlying these reductions in 
efficiency were shown to be comprised of generally weaker 
connections, that is, fewer, smaller or less well-myelinated fibers. 
In the case of nodal global efficiency the paths involved were also 
shown to be on average physically longer and more indirect, that 
is, to pass through more intermediary nodes. In the case of nodal 
local efficiency the paths involved were comprised of weaker, but 
generally physically shorter and more direct paths in the HR^^^ 
infants. But as the degree to which a node is a neighbor of 
another is determined by the strength of the direct connection 
between them, and the nodal global connection strength is 
decreased in the HR'^^^ infants, the physically shorter and more 
direct paths in the HR'^°^ infants represents reduced spatial extent 
of the nodal local networks. Thus, the reduction in nodal local 
efficiency in the HR'^°^ infants also reflects a reduction in long- 
distance connectivity. The reductions in nodal local and nodal 
global efficiency in the HR'^^^ infants are thus not simply indices of 
under-connectivity, but also reflect a more random configuration 
with more limited local connectivity and less direct connections to 
other brain regions. 

Networks with a high degree of spatially local connectivity, but 
with few or no long-range connections, that is, shortcuts, have 
high local efficiency and low global efficiency; networks with a 
high degree of long-range connectivity, but which lack spatially 
local clustering, have high global efficiency and low local 
efficiency. Biological systems in general, and neural networks in 
particular, reliably balance global efficiency with local efficiency, 
having strong local clustering mixed with sufficient long-range 
connectivity to allow rapid communication between distant 
nodes; these have been dubbed 'small-world' properties.^'^'^^'^'' 
The results here suggest that cortico-cortical organization is 
deficient in infants with significant symptoms of ASD in terms of 
both the long-range connectivity that provides for rapid integra- 
tion of information between different brain regions, and the 
spatial clustering that provides for efficient processing for 
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subtasks. Notably, these abnormalities are predominately poster- 
ior, and in regions involved in early processing of auditory and 
visual stimuli, audio-visual integration, language and nonlinguistic 
social stimuli. 

The economical 'small-world' organization of the typical human 
brain is refined over development.^^'^^'^^ The results here might 
represent either delayed or aberrant development of this sort of 
economical organization. Further research is needed to determine 
the developmental nature of the network inefficiencies in infants 
who go on to a diagnosis of ASD and how these inefficiencies 
might relate to brain overgrowth and emerging symptoms.^^'^^ 
The key to understanding development is development itself.^^ 
Longitudinal data spanning infancy and toddlerhood will be 
essential in addressing these questions, as will four-dimensional 
multimodal methods to accurately estimate and analyze changes 
in connection lengths, connection strengths, network organization 
and behavior. 
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